makephi <- function(coords, scale= 10){
	-log(0.05)/(max(dist(coords))/scale)
}

